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Abstract 

In the paper, the principal aspects of the mathematical theory of equilibrium 
thermodynamics are distinguished. It is proved that the points of degeneration 
of a Bose gas of fractal dimension in the momentum space coincide with critical 
points or real gases, whereas the jumps of critical indices and the Maxwell rule are 
related to the tunnel generalization of thermodynamics. Semiclassical methods are 
considered for the tunnel generalization of thermodynamics and also for the second 
and ultrasecond quantization (operators of creation and annihilation of pairs). To 
5jq[ every pure gas there corresponds a new critical point of the limit negative pressure 

below which the liquid passes to a dispersed state (a foam). Relations for critical 
points of a homogeneous mixture of pure gases are given in dependence on the 
concentration of gases. 

"Leontii Sergeevich wrote a play himself!" 
"What for?" asked Natalya Ivanovna with apprehension, 
"Are we running out of plays? Such good plays around. 
And so many. . . " 
"Leontii wrote a modern play!" 
Here the old lady became quite alarmed. 
"We do not start mutinies against the authorities!" 

M. Bulgakov "Theatrical Novel" 



1 Introduction 

In his famous book on thermodynamics [T], Leontovich refers to thermodynamics as a 
phenomenological science, i.e., based on a series of "statements established experimen- 
tally". Beginning with the creation of satellites and experiments in the absence of the 
gravity field of Earth, a new period in physical experimental investigations of thermo- 
dynamical phenomena occurred. Although these experiments are carried out by robots, 
and they cannot react like a human being to any situation and to change the parameters 
of the experiment like an animate observer, this must nevertheless unavoidably lead to a 
reconsideration of the original conception of the great physicists. 

On the other hand, computer-aided experiments have been developed so greatly that 
a new science arose, the so-called molecular dynamics. 
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The problem of the mathematical axiomatization of thermodynamics and probability 
theory was posed already by Hilbert. 

Already Poincare gave the standard definition of probability as the ratio of the num- 
ber of cases favorable for the event to the total number of possible eventd3 and gave 
a counterexample to this definition of probability. This definition must be completed, 
writes Poincare, by the sentence "under the assumption that these cases are equiprob- 
able" ([2], Russian p. 116) and notes that we have completed a vicious circle place by 
defining probability via probability!! 

After this, Poincare writes: "The conclusion which seems to follow from filial is that 
the calculus of probabilities is a useless science, that the obscure instinct which we call 
common sense, and to which we appeal for the legit imisat ion of our conventions, must be 
distrusted." ([2j, Russian p. 116)0 

Thus, first of all, the problem is to define what cases can be most naturally assumed to 
be equiprobable. "We are to look for a mathematical thought," writes Poincare, "where 
is remains pure, i.e., in arithmetic" ([21], Russian p. 14) H 

The relationship between the quantum statistics of Bose gas with number theory 
( "arithmetics" ) was traced by Vershik, one of the greatest contemporary mathematicians, 
who followed Temperley, a famous physicist [3]. Meanwhile, Poincare notes that one still 
cannot reject probability theory, because probability theory is related to thermodynamics, 
and thermodynamics to experiment. 

The author had to study the modern experiments in the most thorough way and prof- 
ited from consultations with the great Russian virtuosi in experimental physics, including 
V. V. Brazhkin, K. I. Shmulovich, A. A. Vasserman, V. I. Nedostup, V. G. Baidakov, 
A. E. Gekhman, and also V. A. Istomin, the greatest expert in applied thermodynam- 
ics and the theory of gas mixtures. It was necessary to study the works of Profes- 
sor D. I. Ivanov especially thoroughly; Professor Ivanov discovered a discrepancy between 
the Widom-Kadanoff- Wilson theoretical critical exponents ([1], [5], [6]) with the experi- 
mental ones ([7], [8]) and carried out classical measurements in thermodynamics. 

Birkhoff, a great mathematician, devoted the first chapter of his book "Hydrodynam- 
ics" to paradoxes in hydrodynamics. In the strict mathematical understanding, a paradox 
is a counterexample to a system of axioms. 

First of all, the system of axioms we have constructed includes a modified analog of 
the "equidistribution law" [9]. 

Our main axiom is the second one, which claims the (almost) uniform distribution of 
densities in a vessel. This axiom leads to number theory, where, instead of presenting 
all the theorems, we address the intuition of the physicists who are acquainted with the 
Bose-Einstein distribution (this distribution is closely related to number theory). This 
is sufficient for a new conception of the notion of ideal gas. Further we use the third 
axiom, claiming the existence of the Zeno-line for pure gases, which enables us to pass 
from ideal gases to unperfect ones. To explain the Maxwell rule and the jumps of critical 
indices at a critical point, as compared with the classical ones, we use the axiom of tunnel 

1 "La definition, dira-t-on, est bien simple: la probability d'un evenement est le rapport du nombre de 
cas favorables a cet evenement au nombre total des cas possibles." 

2 "On est done reduit a completer cettc definition en disant: '. . . au nombre des cas possibles, pourvu 
que ces cas soient egalement probables.' Nous voila done reduits a definir lc probable par le probable." 

3 Poincare presents a series of contradictions in probability theory, including Bertrand's paradox. 

4 "La conclusion qui semble resulter de tout cela, e'est que le calcul des probabilites est une science 
vaine, qu'il faut se defier de cet instinct obscur que nous nommions bon sens et auquel nous demandions 
de legitimer nos conventions." 

5 "II nous faut chercher la pensee mathematique la oil elle est restee pure, e'est-a-dire en arithmetique." 



quantization of the classical equation of state, which thus encloses the system of four 
axioms. 

The theory of ideal liquid treated as an incompressible one is a certain new model in 
equilibrium thermodynamics. 

2 Heuristic considerations. 
The role of small viscosity 

Following Clausius, experts in molecular physics usually argue by proceeding from the 
symmetry of the motion of a molecule averaged in all six directions. In the scattering 
problem, we use the principle of symmetry in all directions, which is standard in molec- 
ular physics, but apply it to define not the mean free path, but other molecular physics 
quantities. Therefore, the fraction of all particles that moves head-on is 1/12. There are 
three such directions; hence, one quarter of all molecules collided 

For the interaction potential, we consider the Lennard- Jones potential 

where e is the energy of the depth of the well and a is the effective radius. 

In the absence of an external potential, the two-particle problem reduces to the one- 
dimensional radial-symmetric one. As is well known [10J, two quantities: the energy E 
and the momentum M are conserved in this problem. In the scattering problem, it is 
convenient to consider, instead of the momentum M, another preserved constant, namely, 
b = y/mp, where p is the impact parameter, in such a way that 

M = VEb. (2) 

By solving with respect to energy E the well-known relation 



p 2 M 2 

E = P - + — +u{r) ) (3) 
m 2mr z 



we obtain the attracting Hamiltonian H: 

H= r ?l(2m)+u(r) a<r<b 
1 — cr /r z 

The repulsing Hamiltonian is separated from if by a barrier. Repulsive particles put 
obstacles in the way of particles of the Hamiltonian H, by creating "viscosity" . 

As the temperature decreases, the barrier height increases up to the value E CT = 0.286£, 
and then starts decreasing (see Fig. 1). According to rough energy estimates [11] , for lesser 
temperatures, an additional barrier must arise when clusters form. This barrier can be 
given for carbon dioxide gas by micelles and for neutral gases and methane by germs of 
droplets, i.e., three-dimensional clusters that contain at least one molecule surrounded by 
other molecules (a prototype of a droplet) 0. 



6 The arguments put forward Clausius concerning symmetry which were applied by Clausius to evaluate 
the free path length and repeated here by the author are quite approximate. However, these arguments 
do not influence on the values of ratios of the form Tb/T ci . 

7 By a "barrier" we mean an obstacle for the collision of particles, namely, an "envelope" of surrounding 
particles that defends the given particle from an immediate blow. In mathematics, a domain is an open 
region containing at least one point. 
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Figure 1: The values of E(b, r) for diverse values of the target parameter. 



However, to study the penetration through the barrier of the incident particle, we must 
plot E along the y axis and turn the wells upside down. Then the minimum becomes the 
barrier and the maximum becomes the depth of the well. 



Figure 2: The trap for a fictitious particle in the the center of mass coordinate system. 
The radius vector r of the /z-point is marked on the abscissa axis. The particle falls from 
the left from the point r = p, where p stands for the target parameter. In the original 
problem ([T]) (before the change of variables M = \fEb) the relation r = p is attained as 
r — > oo. 

A dimer can be formed in a classical domain if the scattering pair has an energy equal 
to the barrier height, slipping into the dip in "infinite" time and getting stuck in it as a 
result of viscosity (and hence of some small energy loss), because this pair of particles, 
having lost energy, hits the barrier on the return path. If the pair of particles has passed 
above this point, then the viscosity may be insufficient for the pair to become stuck: 
such a pair returns above the barrier after reflection. Therefore, only the existence of a 
point E = -E max plus an infinitesimal quantity, where E max is the upper barrier point, is 
a necessary condition for the pair to be stuck inside the dip; -E^ax * s the height of the 
maximum barrier. 
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We can compare the values T cr with the values E™ ax in the table below. 



Substance 


e,K 


T cr /4 


E cr ■ e/ k 


Ne 


36.3 


11 


10.5 


Ar 


119.3 


37 


35 


Kr 


171 


52 


50 


N 2 


95.9 


31 


28 


CH A 


148.2 


47 


43 



Above the value E B = 0.8e, the trap disappears. At the value 0.286e, the depth of 
the trap is maximum and corresponds to T cr = ^rp. For neon and krypton, as can be 
seen from the table, the concurrence is sufficiently good. Because T B = 3.2e/k, it follows 
that Tb/T ct = 2.7, which corresponds to the known relation of "the law of corresponding 
states" [12] . 

The temperature corresponding to AE^/k, is the temperature above which dimers do 
not appear. Exactly this is what we call the Boyle temperature (in contrast to [S]). 

In fact, an application of the Clausius approach to the pairwise interaction gives a 
pairwise interaction with respect to the Lennard- Jones potential for two Gibbs ensembles 
of noninteracting molecules. This leads to the presence of a small friction for a single pair. 

The difference E max — E min is equal to the energy needed for a particle lying at the 
bottom of the potential well to overcome the barrier. The value E max corresponds to 
the temperature given by E max = RT, where R stands for the universal gas constant. 
According to graph 2, E min corresponds to the energy PV. Therefore, E min / E max < 1 is 
the compressibility factor, Z = PV/RT. The temperature at the point E m[n = E mSuX is 
equal to the Boyle temperature. 

The dressed or "thermal" potential (p(r) is attractive [13] . In addition, because the vol- 
ume V is a large parameter, it follows that, if the quantity <p(r) = (Nea/\^V)U(r/\^V), 
where U(r/y/V) is a smooth function and N stands for the number of particles, is ex- 
panded in terms of 1/ \/V , then 



Expanding 



(r - r ) 2 (r + r ) 



2 



Cl + r*=^-^- + ^f^, (6) 

where C\ = r%, we can, just as in [13], separate the variables in the two-particle problem 
and obtain the scattering problem for pairs of particles and the problem of their cooper- 
ative motion for r% + r%- The term CirjyfV does not depend on this problem and the 
correction [aej vV)NO(l/V) is small. 

Then, in the scattering problem, an attractive quadratic potential (inverted parabola 
multiplied by the density or, to be more precise, by the concentration, which we denote 
by the symbol p as well, because the target parameter does not occur below) is added to 
the Lennard- Jones interaction potential. 

For this problem, we can find just as in ©-fll]), for all p = N/V, a point corre- 
sponding to the temperature at which the well capturing the dimers vanishes, and thus 
determine the so-called Zeno-line. It is actually a straight line (up to 2%), on which 
Z = E miQ / E m3iX = 1 (i.e., an ideal curve). 

Let us clarify this fact in more detail. 



We can treat the repulsing potential as a potential creating a small viscosity. 
Let us find the total energy of the attracting Hamiltonian, 



E 



rav 



2(1 -b 2 /r 2 



1 -b 2 /r 2 ' 



$(V) 



u[r) — pr . 



The first term is negative for r < b and the other term is positive for b > r > a (i.e., the 
more is the speed, the less is energy). The mean speed is temperature. 
Let us make the change of variables 



r 
a 



b 
a 



and get rid of a. In what follows, we omit both the tilde and the prime. 

For a given b, the minimum n and the maximum r2 (see the graph no. 1 in [14J) are 
defined by the relation 

(T) 



dr 

This gives E max and E min . These values coincide at some point b = b , and hence 

d 2 E 



dr 2 







(8) 



at the point r , i.e., -E max = E min , and this is the very Zeno-line. 

Let us construct the curve Z min = E min / E max minimal with respect to the target 
parameter as a function of p. Let us find the point Z = E min / -E max for -E max = E^ ax and 
find the corresponding point on the curve Z min (p). This point is equal to Z cr = 0.29, i.e., 
to the critical value of the compressibility factor Z for argon. 



Zeno-line 




Figure 3: The binodal, the Zeno-line, and the curve Z m i n . This heuristic binodal does 
not coincide with the experimental one, whereas the Zeno-line and Z CT are close to the 
corresponding experimental curves. 

In order to obtain a binodal according to some "heuristic principle," we must subtract 
the curve Z min (p) [15] from the Zeno-line. This gives the graph shown in Fig 3. 

3 Quantum case 

In an example describing the creation of a dimer, it is shown that, for T = T c , one degree 
of freedom becomes " frozen" , and we obtain two degrees of freedom rather than three. 



For a dimer with T > T c , if the oscillational degrees of freedom are taken into account, 
then the number of degrees of freedom becomes equal to 6. Two degrees of freedom are 
obtained under the assumption that the oscillational degrees of freedom of the dimer are 
also "frozen" at T = T c . If we suppose this heuristic supposition for the quantum case, 
then, for T < T d , both dimers with two degrees of freedom and dimers with six degrees 
of freedom are created. This corresponds to the two-liquid Thiess-Landau model [67]. In 
this case, the dimers with two degrees of freedom give the A-point and the dimers with 
six degrees of freedom give superfluidity. Indeed, in the two-dimensional case we have 

^2T l°° Jd^ + T_ I™ e^jdj | Q (T-T d 



T d J e* - 1 T d J (ef-M/r _ 1)2 ^ Td 

and we obtain a logarithmic divergence at the point £ = for fi — > 0. 

The case in which N is not so large as it is in statistical physics, i.e., the so-called 
mesoscopic state (see [68]), can also be of interest for us. In this case, let us use Fock's 
idea for the Hartree equations, which lead to the Hartree-Fock equations, which agree 
well with experiments for a not too large number of electrons in atoms. 

Namely, we consider the single-particle equation of the mean field (a self-consistent 
field) and apply (to the resulting "dressed" potential) the procedure of transition to the N- 
partial Schrodinger equation with a dressed potential, just as we proceeded above for the 
operator j-A. Here we can consider two ways of investigation. The first way is the way 
used by Fock and which leads in the semiclassical limit to the Thomas-Fermi equations 
for the dressed potential. Another way is to consider the Hartree temperature equations 
(see [69]) and to obtain the Thomas-Fermi temperature equations in the classical limit. 

Since the quantity T d is small, it is easier to use the first way and to find the "dressed" 
potential. 

Let V(q — q') be a pairwise interaction potential such that j\V(r)\dr < oo. The 
dressed potential W(q) is given by the formula 

W(q) = U(q) + J V(q-q')\ij(q')\ 2 dq', 

where U(q) stands for the external potential and ip(q') for an eigenfunction of the 
Schrodinger equation which depends on the "dressed" potential and is thus an equation 
with a "unitary" nonlinearity. The expansion of the equation in powers of h can be found 
by the method of complex germ up to 0(h k ), where k is an an arbitrarily large number^ 
(see [69], where system (63) defines a complex germ; see also [40], [75-80]). 

Further, we consider two asymptotic formulas for the number of eigenvalues of the 
one-particle Schrodinger operator. 

The 2.D-dimensional phase space is partitioned into a lattice, and the number Gj is 
determined by the formula 

G, = (9) 



This relation was obtained by Weyl. 

In the mathematical literature, the same notion for the number of eigenvalues of 
the Laplace operator with the Dirichlet conditions on the boundary of a D-dimensional 



8 For U(q) = 0, one obtains Bogolyubov's famous equation [70]. The creation of dimers leads to 
the ultrasecond quantization, i.e., to the operators of creation and annihilation of pairs. This makes it 
possible to satisfy the boundary conditions in a capillary (46-48, 71-73). 
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volume, which are less or equal to a given A, is called the Courant formula: 



P(A) 



y m D/2 X D/2 



(1+0(1)) 



as A — > oo 



(10) 



T(D/2 + l)(27r) D / 2 h D 



where D stands for the dimension of the space, because the spectral density has the 
asymptotic behavior. 
The asymptotics 



is a natural generalization of formula (TTTJT) . 

To reconcile the notion of Bose statistics which is given in [9] with symmetric solutions 
of the iV-particle Schrodinger equation, i.e., of the direct sum of iV noninteracting Hamil- 
tonians corresponding to the Schrodinger equation, and the symmetric solutions of their 
spectrum, it is more appropriate to assign to the cells the multiplicities of the spectrum 
of the Schrodinger equation in the way described in [81]. 

Consider the nonrelativistic case in which the Hamiltonian H is equal to p 2 /(2m), 
where p stands for the momentum. 

The comparison of Gi with the multiplicities of the spectrum of the Schrodinger equa- 
tion gives a correspondence between the eigenfunctions of the iV-partial Schrodinger equa- 
tion that are symmetric with respect to the permutations of particles and the combina- 
torial calculations of the Bose statistics that are presented in [9]. 

A single-particle ^-function satisfies the free Schrodinger equation with the Dirichlet 
conditions on the vessel walls. 

Using precisely this very correspondence, we establish a relationship between the Bose- 
Einstein combinatorics [9], the definition of the iV-particle Schrodinger equation, and the 
multiplicity of the spectrum of the single-particle Schrodinger equation. 

The spectrum of the single-particle Schrodinger equation, provided that the interaction 
potential is not taken into account, coincides, up to a factor, with the spectrum of the 
Laplace operator. Consider its spectrum for the closed interval, for the square, and for 
the D- dimensional cube with zero boundary conditions. This spectrum obviously consists 
of the sum of one- dimensional spectra. 

On the line we mark the points i = 0,1,2,... and on the coordinate axes x, y of the 
plane we mark the points with x = % = 0, 1, 2, . . . and y = j = 0,1,2, ... . To this set of 
points we assign the points on the line that are positive integers, I = 1,2, ... . 

To every point we assign a pair of points, % and j, by the rule % + j = I. The number 
of these points is ni — I + 1. This is the two-dimensional case. 

Consider the 3-dimensional case. On the axis 2 we set k = 0,1,2, ... , i.e., let 



It can readily be seen, for the D-dimensional case, that the sequence of multiplicities 
for the number of variants i = Ylk=i m k-> where are arbitrary positive integers, is of 
the form 




as j — > 00 




i+j+k=l 




<k{D) 




for D 




(12) 



s 



Nt = N, eJ2 Qi(D)Ni = E. (13) 

i=l i=l 

The following problem in number theory corresponds to the three-dimensional case 
D = 3 (cf. [9]): 

f> = AT, ef;fi^^ = M, M=Z (14) 

i=l i=l 

In the manual [9], following Einstein, a passage to the limit is carried out as iV — > oo, 
which enables one to pass from sums to integrals. Then, in the section "Degenerate Bose 
gas," a point is distinguished which corresponds to the energy equal to zero. This very 
point is the point of Bose condensate on which excessive particles whose number exceeds 
some value N d ^> 1 are accumulated at temperatures below the so-called degeneracy tem- 
perature Td- The theoretical discovery of this point anticipated a number of experiments 
that confirmed this fact not only for liquid helium but also for a series of metals and even 
for hydrogen. 

From a mathematical point of view, distinguishing a point in the integral is an incorrect 
operation if this point does not form a 5 function. 

The author obtained a more precise expression for the Bose-Einstein but for finite 
values N d 3> 1 (see, e.g., [82]). 

Thus, we consider the case in which N ^> 1, but n is not equal to infinity, on the 
physical level of rigor [82]. For the final parastatistics we have 

n .-_^ k + 1 n --^ (15) 

exp{^}-l ex P {(* + l)^}-l' "> Gj' [Lb) 

In our case, we have k = N d , and the point of condensate is e = 0. 

By (jnj), it is clear that Gj is associated with the D- dimensional Lebesgue measure 
and, in the limit with respect to the coordinates Aqj, gives the volume V in the space of 
dimension 3 and the area Q in the space of dimension 2. The passage with respect to the 
momenta Apj is also valid as N — > oo and /i > 5 > 0, where 5 is arbitrarily small. 

Expanding ( II 5p at the point Eo = in the small parameter x = (//iV^/T^, where 
Nd stands for the number of particles corresponding to the degeneration and T d for the 
degeneracy temperature, and writing £ = —fi/T d , we obtain 



N d + 1 \ e^-l-(iV d + i)(e«-i) 



exp{^} - 1 exp{(iV d + 1)=*} - 1 j (e* - l)(e 2 ^ - 1 



JV d l + f + 2j + ff + ... _ Nd (. , 11 



tl± l_l_±± x 2_ 0191:c 3_ (16) 



n 



2 i + £ + s! + 4 + 2 V 3 24 

For example, if x — > 0, then no = N d /2, and hence the number no in the condensate 
at T — Td does not exceed N d /2. If x = 1.57, then no ~ Certainly, this affects the 

degeneracy temperature, because this temperature can be expressed only in terms of the 
number of particles above the condensate, iVd, rather than in terms of the total number 
of particles N d (which is equal to the sum of N d and of the number of particles in the 
condensate). 

Write M = E d /£i, where Si stands for the coefficient in formula (iTTj) for j = 1. Let 
us find E d , 



2m 



de 



E d = I , (17) 



\P\ It 
e -2^l 1 d — 1 



where 

\p\ 2 dp 1 ...dp D dV D 

d£ = 2^ — {2^r — • (18) 

Whence we obtain the coefficient a in the formula, 

E d = aT 2 d + \(l + D/2)Y(l + D/2). (19) 

To begin the summation in ( fi~3l) at the zero index (beginning with the zero energy), 
it is necessary to rewrite the sums (fT3|) in the form 

oo oo 

Y^N l = N, e^2{ qi (D)-l)Ni = E-eN. (20) 

i=0 i=0 

The relationship between the degeneracy temperature and the number of particles 
above the condensate for /i > 5 > (where 5 is arbitrarily small) can be found for D > 2 
in the standard way. 

Thus, we have established a relationship between Gi in formula ([9]) (which is combina- 
torially statistical) and the multiplicity of the spectrum for the single-particle Schrodinger 
equation, i.e., between the statistical [9] and quantum- mechanical definitions of Bose par- 
ticles. 

Consider the two-dimensional case in more detail. There is an Erdos' theorem for a 
system of two Diophantine equations, 

oo oo 

Y^n 1 = n, j2 iN * = M - ( 21 ) 

i=l i=l 

The maximum number of solutions of this system is achieved if the following relation is 

satisfied: 

N d = c-'M 1 / 2 log M d + aM X J 2 + o(M d 1/2 ), c = vr y/2/3, (22) 

and if the coefficient a is defined by the formula c/2 = e _ca//2 . 

The decomposition of M d into one summand gives only one version. The decompo- 
sition M d into M d summands also provides only one version (namely, the sum of ones). 
Therefore, somewhere in the interval must be at least one maximum of the variants. Erdos 
had evaluated it (J22D (see [83]). 

If the number iV increases and M is preserved in the problem (I2ip . then the number 
of solutions decreases. If the sums (1211) are counted from zero rather than from one, i.e., 
if we set 

oo oo 

Y^m = (M-N), ^iV i = iV, (23) 

i=0 i=0 

then the number of solutions does not decrease and remains constant. 

I shall try to explain this effect. The Erdos-Lehner problem [84] is to decompose M d 
into N < N d summands. 

Let us expand the number 5 into two summands. We obtain 3 + 2 = 4 + 1. The total 
number is 2 versions (this problem is known as "partitio numerorum"). If we include 
to the possible summands, we obtain three versions: 5 + = 3 + 2 = 4 + 1. Thus, the 
inclusion of zero makes it possible to say that we expand a number into k < n summands. 
Indeed, the expansion of the number 5 into three summands includes all the previous 
versions, namely, 5 + + 0, 3 + 2 + 0, and 4 + 1 + 0, and adds new versions, which do not 
include zero. 



in 



In this case, the maximum does not change drastically [84]; however, the number of 
versions is not changed, namely, the zeros, i.e., the Bose condensate, make it possible 
that the maximum remains constant, and the entropy never decreases; after reaching the 
maximum, it becomes constant. 

Let us turn to a physical definition. 

Note first that, without changing the accuracy of the quantity whose logarithm is 
evaluated, we can replace logM^ by (1/2) log(Nd/Q)- Then 

^- 2f >«l Q + „f*Y (24) 



c-n g(N d /Q) + a \Q 

In our case, N d /Q corresponds to the number of particles above the condensate. 

According to formula (TO?]) , in the two-dimensional case we must set 7 = and find 
the coefficient a. Then formula f)24p gives us a relationship between and due to the 
fact that the number of particles in the condensate is o(Nd). 

In fact, we have proved that there is a gap between fi > 5 > and \i = 0. 



4 Main axiom and theorem for the classical gas 

Consider the Maxwell distribution 

w{p) = (2 Jt)3/ 2 exp { ~ £r }• (25) 

As is well knowrj^l, it is associated with the potential 

V7r 3/2 T 5/2 m 5/2 ,00 

tt = — — - — — — / t 3/2 e- {t -T } dt. (26) 



1X2 + 1/2) J 

Now consider an ultrarelativistic gas in which the kinetic energy is proportional to \p\. 
In this case, 

»®= (27) 

where c is the velocity of light. For this gas, the potential Q is of the form 

Since = — |p and P = — |p, it follows that, in both the first and the second case, we 
obtain the equation of an ideal gas 

PV = NT. (29) 

In thermodynamics, there is an important notion, namely, I, the number of the degrees 
of freedom of the molecule. The energy of the molecule is 

v 2 



3 We have kept the integral especially without evaluating it. 



1 1 



The parameter I (the number of degrees of freedom) is the coefficient at the potential Q. 
It does not influence the equation of ideal gas ( 1291) . This coefficient is related to the so- 
called "equidistribution law" mentioned above. Possibly, the law is true in some interval 
of temperatures. However, as the temperature reduces, some degrees of freedom can be 
"frozen." It turns out that it is much more convenient to introduce the parameter a 
presented below. 

Assuming that the value of I depends on the momentum, we replace the above relation 
by the following one: 

\v\ 2+a 

£7 = — , (30) 

where po stands for some typical momentum. 

We can say that both the average momentum and the number of the degrees of freedom 
depend on the temperature, and hence the number of the degrees of freedom depends on 
the momentum. We assume that this dependence is the simplest one, namely, it is a 
power dependence (the parameter a of in the exponent characterizes the molecule). Here 
the potential Q takes the form 

-7r 1+ 7 m 2 +7 f°° 



o=^——VT 2+ ~<-^ -/ t^e-^dt, (31) 
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7 = 1 = , 

' 2+a 2+ct 

where A stands for a constant corresponding to a given molecule and depending on its 
mass (however, we try to avoid the mass by passing from the density to the concentration). 
It can readily be seen that, in this case, relation f[2"9j) remains valid for an ideal gas. 

Thus, we replace the parameter of the integer degrees of freedom by some continuous 
parameter that characterizes the given molecule. In principle, this parameter 7 is of the 
same physical origin as the number of the degrees of freedom. But since it is continuous, 
it takes into account more details of the spectrum of the molecule. 

Moreover, below we shall see the following remarkable fact: every molecule is charac- 
terized by its own value of 7 which does not depend on temperature and depends only on 
the dimensionless value Z CT , which is its critical compressibility factor. 

One must not think that this relation (formula fj30|) ) holds for a single molecule. It 
can be treated only as a result of averaging over the Gibbs ensemble of noninteracting 
molecules that corresponds to some temperature (as we shall see below, this ensemble 
corresponds to the critical temperature T cr ). 

Roughly speaking, the number of degrees of freedom is related to the number of 
atoms [9], whereas the number of atoms is related to the dimension of the space. For this 
reason, the reduction of the number of degrees of freedom, as the temperature reduced, 
is equivalent in a Gibbs ensemble of noninteracting molecules to the reduction of the 
mean value of the space dimension. On the other hand, one can say that the reduction 
of the number of degrees of freedom is analogous to the reduction of the dimension of 
the space of oscillators (from three-dimensional to zero-dimensional oscillators), whereas 
the latter reduction is equivalent (see Sec. 5) to the halving of the dimension of the 
momentum space. However, for the value of the potential Q, increasing of the exponent 
of p in energy is equivalent to the reduction of the fractal dimension when integrating with 
respect to p. Therefore, the passage to E{p) = ap a , which is regarded as an equivalent 
notion of fractal dimension, is related to averaging of the Gibbs ensemble of noninteracting 
molecules with respect to the number of active degrees of freedom. Therefore, the formal 
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change of variable in manuals Q60J, p. 212) and handbooks [61] (see the section "Integral 
representations" in [61] and also the book [62]) has a well-founded physical meaning 
of continuous reduction of a fractional value of the number of degrees of freedom of a 
molecule as the temperature reduces (in the mean) for a Gibbs ensemble of noninteracting 
molecules. 

A. V. Chaplik suggested another interpretation of this law. Since it turns out that 
a law of the form ap a enables us to describe the experiment, and since we assume that 
this is an ideal gas of particles with this dependence, E{p) = ap a , it follows that these 
are quasiparticles, i.e., collective excitements (rather than individual molecules). The p- 
momentum of the quasiparticle is preserved, because the interaction between the particles 
is no longer valid. This trick is customary in physics. 

However, in mathematics, the notion of Gibbs ensemble corresponds to the notion of 
the Kolmogorov complexity, and therefore, for a mathematician, it is more natural to 
relate the formula E(p) = ap a to this vary complexity. However, the main point is that 
this generalizes the notion of equidistribution to lower temperatures. 

Now we can rigorously formulate the axiom of thermodynamics corresponding to the 
approximate conservation of the gas density (this corresponds to the physicists' statement 
in equilibrium thermodynamics: "the density is homogeneous in a vessel"; physicists 
consider equilibrium thermodynamics as a separate discipline, after which fluctuation 
theory is also considered separately). 

A physicist thinks that the density of particles inside the vessel is constant up to fluc- 
tuations (of the order of the square root of the number of particles in a small subvolume). 
This means that the physicist counts the number of particles in every small subvolume 
containing approximately one million of particles, and the density does not depend on 
the permutation of the numbers of particles. However, this is arithmetics, in which the 
sum does not depend on the permutation of summands (in contrast to the Boltzmann 
statistics in which a permutation of particles gives a new state). Even if we assume that 
the experimenter had indexed all particles during the previous measurement, this exper- 
imenter cannot control what index corresponds to a given particle at the next moment 
of measurement, and the indexing can be carried out again to compute the density. This 
is a fact in arithmetics, which is the very base of a special branch of science, namely, of 
analytic number theory. 

Main mathematical axiom of thermodynamics. Consider a vessel of volume V 
containing N > 10 19 identical molecules corresponding to the parameter 7 = 70 . Consider 
a small convex volume of size V$ containing Ng particles, where N$ is not less than 10 6 . 
Let P be the probability of the deviation 



for any volume of size V$ inside the vessel. The claim of the axiom is that the probability P 
is sufficiently small. Namely, 



This implies that the numbering of particles in the volume V5 is arbitrary and does 
not affect the density. We can rearrange the numbers, and this also does not affect the 
density, because the rearrangement of summands does not change the su 

10 The set of initial data for an TV-particle dynamical system can be said to be chaotic if a relation 
of the form (|32|) arises at a sufficiently large time interval. The weakening of condition (f32| does not 
influence the validity of relation ([M)) . 



N 5 ± ^N~ 5 




(32) 



Remark 1. Condition (132]) shows that the theory presented below cannot be applied for 
insufficiently many particles. This condition also shows what is the accuracy that can 
be expected from the relations obtained below. The weakening of condition ( 1321) plays 
no fundamental role, and is manifests itself only at the level of the expected accuracy of 
these relations. 

Thus, we have already proved that the approximate conservation of density inside a gas 
volume implies that the particles are independent of the numbering, and hence are Bose 
particles. Therefore, the potential f2 7 ( I3~il must be replaced by the Bose gas potentiate, 

where Li 7 (-) stands for the polylogarithm, Li 7+2 (1) = C(7 + 2), and ( stands for the 
Riemann zeta function. 

Thus, in our case, the Bose gas also depends on another parameter, 7. We refer to 
this parameter as the fractional {"fractal") dimension in the momentum space (cf. [9], 
Sec. Degenerate Bose gas, footnote), the parameter could be called more precisely the 
parameter characterizing the spectrum of the molecule. 

Remark 2. In the statistics of Boltzmann, Shannon, and Kolmogorov, 4 + 1 and 1 + 4 
are different variants. In arithmetic, they constitute one variant. For 7 > 0, the Bose 
statistics is related to number theory. In the physical literature, this was first noted by 
Temperley [3J; in the mathematical literature, this was considered in detail in [T6]-|18]. 

Remark 3. Note that, just as for the standard potential f2 7 (1311) . three-dimensional Bose 
gas corresponds to 7 = 1/2, whereas 7 = 2 corresponds to ultrarelativistic gas. The latter 
is associated with the famous Planck distribution at fi = for the black-body radiation 
(cf. [9], Sec. 60). 

Theorem 1. The point of degeneracy of the Bose distribution at fj, — for molecules of 
characteristic 7 coincides with the critical points of these molecules up to normalization. 
Therefore, the dimensionless quantity 

V P 

v cr- 1 cr 



z 



where R stands for the gas constant, coincides with the dimensionless point of degeneracy 
of an ideal Bose gas of fractal characteristic 7. Hence, 

von,, /on^ _ ((7 + 2) 



Zc *-\T^/^)^ -tWTTy (34) 

which is a consequence of the number-theoretic interpretation of the Bose distribution. 

This is a property of the new ideal gas without interaction. Apfelbaum and Vorob'ev [19 
computed the critical isotherms in the {P, \^}-plane for different gases, where the value 
of 7 is determined by the above formula, and the critical isotherm of the Bose gas with 
the evaluated characteristic 7 was computed by the relevant formulas for an ideal Bose 
gas El (see Figs. 4-5). 



n The Bose gas is usually regarded as a quantum gas; however, it is related to number theory (see 
Remark 2 below) , and number theory is related to the main axiom. Therefore, the potential f2 7 (|33|) can 
be applied to classical gas. In number theory, V = 1. To make the perception by the physicists more 
convenient, we keep the multiplier V (cf. (|4|)- (p8| below). 

12 This gas is more general than that in [20] and corresponds to the fractal dimension 7 in the momentum 
space. In [20 , 7 — 1/2. There is no classical gas of this dimension. 
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Figure 4: (a) Isotherms of pressure for the van der Waals equation are shown by continuous 
lines. The small circles show the corresponding lines computed with 7 = 0.312 for <p(V) = 
V (i.e., ideal Bose gas), Z CT = 3/8. 

(b) Isobars of density for the van der Waals equation are shown by continuous lines. Line 1 
is the binodal curve. The small circles correspond to Bose-Einstein isobars for 7 = 0.312. 



At first glance, it looks as if the notion of new ideal gas leads to an alteration of the 
famous relation 

PV = NT (35) 

(which, moreover, served as an analogy for the main economics law, the Irving Fisher 
formula; which is used to calculate the "turnover rate" of capital [2T]). This would be 
surprising indeed. However, this is not the case. The relation PV = NT or, equivalently, 
PV = RT (because the number of particles in the vessel remains the same) defines an 
imperfect gas and, in contemporary experimental thermodynamical diagrams, it is called 
the Zeno-line or, sometimes, the ideal curve, the Bachinskii parabola, etc. 

On the diagram (p,T) for pure gases, this is the straight line Z — 1. The line is a 
most important characteristic feature for a gas which is imperfect. Since, for imperfect 
gases, it has been calculated experimentally and is an "almost straight" line on the (p, T) 
diagram, it follows that the Zeno-line is determined by the two points Tb and ps called 
the "Boyle temperature" and the "Boyle density." In contrast to Z CI , these points are 
related to the interaction and scattering of a pair of gas particles with the interaction 
potential peculiar to this gas, as was shown in Sec. 1 and in other papers of the author 
(see, e.g., [22] )• Therefore, the Zeno-line on which the relations 

PV = NT, -L + — = 1, - = p, (36) 
Pb T B V 

are satisfied, where p is the density (the concentration), is a consequence of pairwise 
interaction, and thus is a relation for an imperfect gas. 

The correction related to the existence of the Zeno-line leads to a differential equa- 
tion |23J whose numerical solution yields an alteration to the gas spinodal for every par- 
ticular pure gas. For argon and CO2, this modification is shown in Fig. ^' 



13 By the heuristic considerations related to the scattering problem (Sec. 1), the final point [23] of the 
gas spinodal is equal to Z = 3/2 Z CI , and the spinodal can be approximated by a line segment. (Ideally, 
in infinite time, a fictitious particle (a pair) falls to the bottom due to friction, i.e., the orbit of this 
particle is circular, and thus one degree of freedom disappears. This means that the compressibility 
factor Z = 0.444 at the point E™ ax is reduced by the factor 2/3, i.e., Z CI = 0.296.). This makes it 
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Figure 5: (a) Isotherms for argon. The continuous lines correspond to the experimental 
data, and the line formed by small circles is constructed according the isotherm of ideal 
Bose gas. 

(b) Isotherms for water. Z cr = 0.23. 

(c) Critical isotherms for copper. Z CT = 0.39. 



The distribution of number theory, as opposed to the Bose-Einstein distribution 
does not contain the volume V. Let us consider the distribution of number theory multi- 
plied by unknown function ip l0 (V) which does not vary for 7 > 7o and T < T cr . Then it 
follows from (|36|) that 



r(7o + 2) J e- K e £ -l' ^ dV ' 

V4 (V) Li 70+2 (y ) = ^ 7q+ — ^— +1 , p = ^, y = e K , « = (37) 

The differential equation for ip l0 is 

^(V0Li w+a (y)_ 1 r = : R_ (3g) 



f 70 (V) Li 70+1 (y) ' p 

See Fig. 6. 



possible to construct two points near P = 0, Z = Z CI by the theories of a new ideal (Bose) gas and by 
the fact that the chemical potential of the gas is equal to the chemical potential of an ideal liquid, and 
thus to approximately reconstruct the Zeno-line. 
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Figure 6: The dotted line shows the Zeno-line Z — 1. The bold line is the critical 
isotherm of an imperfect gas (argon) calculated theoretically; the thin lines correspond to 
the isochores of an imperfect gas for T < T cr . Their initial points lie on the gas spinodal. 

5 Ideal liquid 

Let us now pass to the notion of ideal liquid. For an expert in mathematical physics, 
an ideal liquid is an incompressible liquid. In our mathematical conception of thermo- 
dynamics, we shall abide by this definition. In this case, on the Zeno-line on the plane 
{P,Z}, the point P(T,p) in ([36]) is uniquely defined for Z — 1. The isotherm T = const 
is a straight line. The second point is a point of the spinodal. 

As is well known, the passage from the gaseous state to the liquid one is accompanied 
by the entropy drop. Naturally, the entropy, which determines the measure of chaotic 
behavior, is less for the liquid state than for the gaseous state. At the same time, the 
general property of "choosing" a subsystem with the greatest chaoticity among all possible 
subsystems leads to the property of constant entropy of the liquid, which was noted 
both theoretically and experimentally, even if the temperature tends to the absolute 
zero ([63], [&3]) (the entropy tends to log 2). 

In our model of ideal liquid as an incompressible liquid, we suppose in addition that 
the maximum of the entropy on a given isotherm (i.e., as \i — > 0) does not vary when the 
temperature varies (see below Section 6). 

The big thermal potential is of the form 

7T^+%T 1 f°° t 1+ ^dt _ 7r 7+l jL 2(l+ 7 ) T 2+ 7 

" ~ ~ PV " ~ ~ A?m r(2 + 7 ) Jo (e7*)-l ~ 2+l[Z) - { } 

The entropy becomes 

/f)0\ r 2(1+7)7.1+7 1-2(1+7)7-1+7 ,, 

S = -{df) v ^ = ^ + 7) A2(1+7) Li 2+7 (.) Li 1+7 (z)| (40) 

= A 2 ( i +7 ) [( 2 + 7) Li2+ 7 (^) - Li 1+7 (*)-] . 



The maximum at /i = is 

" = ' \ 



^=o=(^) 7+1 (2 + 7 )C(7 + 2)T^. (41) 
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We are interested in 7 < 0. 

Thus, we have two unknown constants, namely, A and the value of the entropy S^o = 
const. We define these two constants from the experimental value of the critical point of 
the liquid phase at the negative pressure (see Sec. 6 below), namely, from the minimum 
point of the pressure for a given simple liquid of the value 7 at this point and from the 
temperature. This point is absent in the van der Waals model. This point is present in 
our model of liquid phase. 

For example, for water, we obtain S^=q = 3.495 and A = 3.74. However, the compu- 
tation is carried out under the assumption that the Zeno-line is a line segment, whereas 
this segment becomes curvilinear for water at low temperatures. 

According to the van der Waals conception, we normalize as follows: 




0.1 j — , — , — , — , — , — , — . — , — , — , — . — , — , — , 
1.0 2.0 3.0 4.0 5.0 6.0 p, 



Figure 7: Experimental graph for the different gases, including for methane, ethylene, 
ethane, propane, n-butane, isopentane, n-heptane, nitrogen, carbon dioxide, and water. 
Each gas is equipped with a particular symbol on the graph. 

Denote by p the solution of the equation 

J-B PB 

where Tb stands for the Boyle temperature and 



n cd = ^ (44) 



- cr 



is a dimensionless quantity, T red = Then 



1 " ^id> Ps d = -■ (45) 



P0_ = l L_ „red _ PB 

B PO 

Hence, the locus of the points of the quasi-ideal spinodals^l is given by the formula 



,r,d = 7i r^^-^zd), (46) 



That is, of the endpoints of the metastable state of the liquid phase. 
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Figure 8: The continuous lines are the experimental isotherms for T > T cr for methane, 
the dotted lines are experimental isochors. The theoretical critical isotherm coincides 
with the experimental isotherm up to Z = 0.29 and is continued by a straight line up to 
the point P — 1, Z — 1/po — 0.14 (see (27)-(29)). Further, at an acute angle, a tangent 
to the experimental isotherm at a point of Zeno-line is drawn. The straight line from the 
point P = 1, Z = 0.14 to the point of tangency is the critical isotherm of the ideal liquid 
phase. It can be seen by comparing the figure with Fig 7, the theoretical isotherm thus 
obtained corresponds to the isotherm of the "law of corresponding states" for the gases 
indicated in Fig. 7. 



where ps and Tb is the Boyle density and the Boyle temperature, respectively. 
Here, for 7 > 0, 

(Trcd) C(2 + 7(T- d )) 
C(l + 7(^ rcd ))' 

Let us recall that 7(T) can be calculated from the algebraic relation 5*^=0 = const. 

If we use the Maxwell condition, see Sec 5 (if Sec 6 is taken into account), which states 
that the transition from gas to liquid occurs for the same chemical potential, pressure, and 
temperature, then we can construct the so-called binodal. The binodal thus constructed 
coincides with the experimental one, in contrast to the van der Waals binodal (Fig. 4 (b)) 
and to the binodal presented in Fig. 3. 



6 Wiener quantization of thermodynamics 
and a jump of critical exponents 

The problem involving the so-called Maxwell rule for the transition "gas-liquid," which 
is a natural complement of the new concept (of phenomenological thermodynamics) con- 
structed above, is solved by using the tunnel (or Wiener) quantization introduced by the 
author already in his 1994 works; see [26] and (27] and also [28], [29], and [30]. We repeat 
here this quantization at a heuristic level. 

We can say that the quantization of thermodynamics is simply called for. Indeed, we 
have the phase spaced in which the momenta are the extensive quantities V and — S, 

15 "A symplectic structure" . 
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The change of variables 
leads to the equation 



and the corresponding coordinates are P and T. The usual quantization is of the form 

Just as in [21], let us invoke an analogy between the Schrodinger equation and the heat 
equation. 

A. The Schrodinger equation corresponding to a noninteracting particle without an 
external field, 

_ ifS$ = (iW) V- (49) 
at 

vb = e^ 

— + (V§) 2 + ihAE = 0. 

In this case, the quantization of the classical Hamilton-Jacobi equation consists in the 
addition of the term ihAE>. 

B. The heat equation 

~ V ~dt = (Z/V)V (50) 
v is the kinematic viscosity. The change of variables 

u = e~~ 

leads to the equation 

^ + (V§) 2 + uAS = 0. (51) 
ot 

The derivatives of this equation with respect to the coordinates are called the Burgers 
equations. In this case, the Wiener quantization consists in the addition of the viscous 
term. 

Remark 4. In this special case, the Wiener quantization coincides with the Euclidean 
quantization well known in field theory. In the general case, it corresponds to the passage 
from the Feynman path integral to the Wiener path integral and is in essence presented 
in detail for physicists in the book [32J. 

In the Burgers equation, for p = dS/dx, a shock wave occurs as v — > 0, i.e., a 
discontinuity of the ^-function type, whereas, in thermodynamics, we have a jump of 
the 6*-function type for the transition "gas-liquid." For the Burgers equation, the rule 
of "equal areas" arises. For the "gas-liquid" transition, the Maxwell rule of equal areas 
arises. In the heat equation, the Wiener quantization of energy yields vD t , where D t 
stands for the Heaviside operator, Dt — Jr. In thermodynamics, the thermodynamical 
potential, the Gibbs energy, is equal to fiN, where \i stands for the chemical potential 
and iV for the conjugate extensive quantity, the number of particles. Hence, 

Of! 

and the role of tinted is played by log — /i, because, under this quantization, the operator 
vii{djdii) corresponds to the Gibbs energy. 



16 Cf. Matsubara Green function, where the role of the imaginary time is played by the parameter 
P = l/T [55]. 
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For us, the one-dimensional case p\ — V and q% = P is of importance. In the general 
case, the focal point (the point of inflection [35J) is of the form q ~ p 3 , i.e., 

P CI ~ (V - V cv ) 3 = ( ??—£r)\ (52) 

V PPcr J 

which corresponds to the classical critical index (the exponent) equal to three. The 
asymptotic solution (as v — > 0) of fl5Tl) corresponding to this point is expressed by the 
Weber function. The critical point itself corresponds to the point x = 0. 
The Weber function is of the form 

-r pea c £ 4 

e'—dC (53) 



'v JO 

If we set x = 0, then the change of variables 

£ = vW /4 

gives a singularity as v — \ of the order of v~ l / A . What does it mean, from the classical 
point of view and classical measurement, for the case in which the condition referred to 
in [20] as the "semiclassical condition" holds (this means that we are outside the focal 
point)? For the Laplace transform, this means that we are in a domain in which the 
Laplace asymptotic method is applicable, i.e., in a domain for which 

*(*) = 4= r d ^ lim < °°> S {i) (P)\v=o + 0. (54) 

The solution p v {x) of the Burgers equation can be evaluated by the formula 

r^EM^ ■ (55) 

After the substitution 

y, 



as x y we obtain 

p u (x) — > x ^o \fu ■ const. (56) 

In our case, the momentum p u (x) is the volume V. 
If the solution of the relation 

x = dS/dp 

is nondegenerate, i.e., 

dp 2 

at the point 

dS _ 
dp 

then in this case, the reduced integral ([55]) is bounded as v — > 0. For this integral to 
have a singularity of order i/ 1 / 4 , we must apply to this integral the fractional derivative 
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Z} -1 / 4 with respect to x. The value of D~ l / A at the function equal to one, .D -1 / 4 !, gives 
approximately x 1//4 . 

In our case, the pressure P plays the role of x, and the volume V plays the role of 
momentum p. Therefore, V ~ P 1/4 , i.e., 



Following M. Green [36], D. Yu. Ivanov, a deep experimenter, poses the following 

question: Why are the deviations from the classical theory in the critical opalescence 
observed within the limits of hundredths of a degree from the critical points, whereas the 
deviations in thermodynamical properties show [7] a nonclassical behavior at a much larger 
distance from the critical point?. Professor Ivanov claims that rather many questions of 
this kind have been accumulated (see, for example, [8]), and all these questions mainly 
deal with the behavior of practical systems. The point is that, from the point of view of 
the developed theory of critical indices [6] , there must be a drastic passage to the classical 
indices outside a neighborhood of the critical point. 

To make Ivanov's question understandable for persons who are not experts in critical 
points, we paraphrase the question for the case of geometric optics, when the sun rays 
are collected by a magnifying glass to a focus. If we had created the special construction 
for the vicinity of the focus in which the paper smoulders, then the experimenter could 
ask why the experiment gives a smooth picture of transition in the double logarithmic 
coordinates and the indices are preserved far away from the smouldering small vicinity of 
the focus. In the present case, the smouldering paper can be compared with the small area 
of opalescence (drastic fluctuations near the critical point) for which a separate theory 
was constructed. At the same time, the special function defining the point within the 
wave theory (like the Weber function) can be continued quite smoothly to a much wider 
domain in which the paper does not smoulder. In the opinion of Ivanov, this fact is much 
more important than the fact that Wilson's theory gives the index 4.82 rather than 4.3, 
which is given by modern experiments. 

Can the experimental index 4.3 be explained in principle in the framework of the 
conception presented by the author? 

Whereas in classical mechanics there is no dependence on the Planck constant h, in 
classical thermodynamics there is a slow dependence of the viscosity on T and p, and thus 
vice versa as well. In our picture, the "stretching" of P cr and T cr in the experiment for 
real gases (Figs. 7 and 8) is greater than in the van der Waals model, which enables us 
to introduce in the "stretching" the parameter v e (i.e., P = v £ V 3 ) and obtain the index 
5 = 4.3 for e = 0.07. Thus, in principle, the answer can be "yes." 

It follows from what was said above that new critical indices arise only due to quanti- 
zation of the conjugate pairs {P, V} and {T, S}. Thus, the relationships between intensive 
quantities can be taken classical, because everything is carried out under the assumption 
of infinitely small viscosity. In this case, one can also pass to another coordinates, to the 
pressure and density. 
Write, as usual, 



Per ~{V~ V cr ) 4 ~ ( 



P ~ Per 




PPcr 




P-P 



cr 







T-T ( 



cr 



P - Per 



P = 



P 



V = 



cr 



cr 



Pcr 



In the classical case, we have Q37J, p. 344) 



P 




v „ e 1/2 , e ~ v 2 . 



In the classical case, we obtain 

p ~ 

In the tunnel quantum case, we obtain j3 = 0.375 (cf. [37], p. 356) in the limit as 
v — y 0, and thus not precisely (the stretching is not taken into account). This is obtained 
for the van der Waals model quantized in the tunnel way. 

An important specific feature of Wiener quantization is that, at the expense of viscos- 
ity, it smoothes the acute angle on the graphs 8 and 9, and thus removes the discrepancy 
between the theory and the experiment. 

The heuristic role of viscosity when penetrating through the barrier is described in 
Sec. 2. As is well known, an exponential penetration through the barrier occurs in the 
standard quantization as well. However, in the Wiener quantization, we restrict ourselves 
to a power-law approximation (see Sec. 6 below). 

If v > 0, then an uncertainty principle arises. Let us proceed with the consideration 
of this principle. 

Thus, in the case of Wiener quantization, the main operators are not the momentum 
operator ih-i^ and the operator of multiplication by x, but the Heaviside operator uD = 
v-i^ and the operator of multiplication by x. The constant v was determined in [38] as 
viscosity. In contrast to the Plank constant, the constant v can smoothly depend on 
the variables of the system. The main distinction of the tunnel quantization from the 
Euclidean and usual quantization is that it is considered up to 0(u k ), where k is any 
fixed number, i.e., it is factorized in 0(v k ). 

But first it is necessary to define the space where these operators act. 

As is well known, the Heaviside operator is related to the two-sided Laplace transform. 
This was shown already by van der Pol and Bremmer in [39]. Introduce a family of 
functions ip(p) to which we shall apply the two-sided Laplace transform, namely, 

POO 

<p(p) = / e-^5(0 
Jo 

and thus these functions themselves are one-sided Laplace transforms, 

POO 

F A S(0= / e" A «H(0^, 
Jo 

for A = p 2 . Denote by the two-sided Laplace transform, 

/oo 
e- x ^(p)dp. 
-oo 

If the functions H(£) are compactly supported and infinitely different iable, then the 
closure of the operator vD with respect to this domain can be carried out in the Bergman 
space. Then the functions *f?(x) = F^tp(p) here become analogous to the ^/-functions 
in the Schrodinger quantization. Moreover, ty*(x) = ^f(x), because these functions are 
real- valued. 

Let us note first of all that the squared function or, equivalently, the squared dispersion 
A f of the operator / is 

(A/) 2 = |(/-/) 2 |, 

and, since / is not selfadjoint, the function (/ — f) 2 need not be positive, and hence 
one must pass to its absolute value. Therefore, the corresponding theorem for generic 
operators fails to hold in general. However, for the operators uD and x on a reduced 
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function space, we obtain |Ai/Z)||Ax| > u/2. It can readily be seen that Weyl's proof 
(which is presented in the comments to §16 of Chap. II in [20]) can easily be transferred 
to the operators vD and x in the above function space. 

Let us note some consequences of tunnel (or Wiener) quantization for the "quantum" 
Bose gas. 

A specific feature of photon gas, which is mentioned in [9], §§ 62,63, is that the number 
of particles in this gas, N, is a variable quantity (rather than a given constant, which is 
the case for ordinary gas). 

Thus, since the number of particles N in thermodynamics is conjugate to the chem- 
ical potential, it follows that, if the number of particles is undefined, then the chemical 
potential can be given precisely, fi = 0, under the assumption that /i and N are tunnel 
quantized and the uncertainty principle holds. 

A contradiction between the conception of the author and the conception of physicists 
going back to Einstein is also removed. In the case of a gas for which iV is fixed, we have 

oo 

^Ni = N (58) 

i=0 

according to the relation in [9], and the chemical potential fi can be a small positive 
quantity. This is obvious, because iVj < N; however, this contradicts Einstein's original 
conception claiming that fi < 0. This contradiction is removed if the relationship of the 
uncertainty principle holds for \x and N, because, if \x = 0, then N can take infinite values 
as well, and therefore the case /i > is impossible. 

Thus, it can be said that both the scaling hypothesis and the hypothesis of tunnel 
quantization do not agree near the critical point with the old thermodynamical conception 
of the four potentials. However, the hypothesis of tunnel quantization does not contradict 
the conception of the four potentials, namely, the hypothesis complements the conception, 
and this works not only near the critical point but also on the entire domain "gas-liquid" 
by agreeing with the Maxwell rule and by removing logical discrepancies in Bose gas 
theory. 

We have explained Wiener quantization, which enabled us to remove some problems. 
Let us now also explain the second quantization. The principal element in Fock's approach 
is the indistinguishability of particles. In our theory, it follows from the original axiom. 
Although there are no natural Hilbert spaces here, in contrast to quantum mechanics, 
we can still obtain correct distinguished representations and limits as h — » (see |40j . 
Chap. 1, Appendix l.A) and then, in view of the indistinguishability of particles, perform 
the second quantization of classical theory by introducing the creation and annihilation 
operators. Certainly, this is possible only under the condition of the principle of indistin- 
guishability of particles in our measurements, which follows from the main axiom, rather 
than from the "identity principle" introduced in [20J. 

In classical mechanics, such operators were introduced in [UJ and (12] on the basis of 
the Schonberg concept (see jl3] and [33]). 

Thus, the contemporary derivation of the Vlasov equation is obtained by applying 
the method of second quantization for classical particles [40J. In this iV — )• oo, 

one obtains a system for which the creation and annihilation operators asymptotically 
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commute, 



,'dU d d\ . 
U ^ q ^=^d- p - p d- q ) U ^ q ^ 

j i j i (i , ^f dV(q,q') d dV(q,q') d\ 
dp dq v(p,q, t) ^ — — + — — Ju(p,q, t)u(p, q, t), 

,'dU d d \ . . . 
v(p, q, t) = -5—5- - P^r )v(p, q, t) (59) 



dq dp dq 

+ y dpd g u(p, g ,t)^ — — + — — )v(p ,q ,t)v(p,q,t), 

where Ufa) is an external field and Vfa,qj) is the pairwise interaction. 

If one replaces u and v by the operators of creation and annihilation u and v in the 
Fock space, then, after this change, system ( 159]) becomes equivalent to the iV-particle 
problem for the Newton system. 

However, according to the mathematical proof, this can happen only for the case in 
which the classical particles are indistinguishable (from the point of view of the main 
axiom) . 

Only in this case the projection from the Fock space to the 3iV-dimensional space of 
N particles gives precisely the system of Newton equations. 
Note that the substitution 



u(p, q, t) = ^^q~Jy^\ v(p, q, t) = y/pfctf)e**<***\ (60) 

reduces system fl59l) to the form 

, (dW l d d \ . 
P^^=[-dq-dp-%)p^^ 

fdW*d d\ f , dV(q,q>)d7i(p',q',t) 

n\P,Q,t) = ^^^r -P^T p(P>9>*) + dpdq — — p(p,q,t); (61) 



dq dp dq J ' J dq' dp' 

we have here 

W\q) = U{q)+ ! dq'V(q,q')p(p',q',t)dp'dq'. 



The first equation of system (45) is the Vlasov equation (see [56J), where p stands for 
the distribution function and W t (q) for the dressed potential (see Sec. 2, formulas ©ff.); 
the other equation is linear, and its meaning is discussed in [59] . 

Note further a tunnel-quantum jump of the index at the points of the spinodal of the 
liquid phase. The classical index of the spinodal is equal to 2, namely, P ~ V 2 , similarly 
to turning points in quantum mechanics. The Airy function corresponds to it. Similarly 
to (IH2]l - (l5i|l . we obtain 

*(*) = 4= He-^dp, lim^<oc, S^(p)\ p=0 ^0. (62) 

The solution p v (x) of the Burgers equation can be evaluated by the formula (1531 . As 
x — > 0, after the change = y, we obtain 

p v {x) — >a;^o • const. (63) 

In our case, the momentum p v {x) is the volume V. Hence, similarly to the consideration 
between the formulas (|55|) and (!56|) . we obtain P ~ V 3 , and the index at the points of 
the spinodal becomes equal to three. 
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Remark 5. It is possible that an experimenter, when considering the approaching of the 
critical isotherm for T > T cr to the critical point p = 0, moves (due to the indeterminacy 
principle) towards increasing values of N, and hence towards increasing density, and 
arrives at the spinodal of the liquid phase. This effect is similar to the accumulation of the 
wave crest which overturns afterwards (a part of the particles outruns the point of creation 
of the shock wave). In this case, the critical index 4.3 passes to the index 3 of the spinodal 
(and this index occasionally coincides with the classical index of the critical point). This 
passage, which is described by the Vlasov equation, was experimentally noticed in [8] and 
in other works. Therefore, the experiments of Ivanov [8] and Wagner ([57]. [58]). where 
the modifications of the critical index S from 4.3 to 3 were obtained when approaching 
the critical point, do not contradict our conception. 

For the creation of dimers, the author of the present paper used the creation and 
annihilation operators for pairs of particles and referred to it as the ultrasecond quanti- 
zation [3n]. Experimenters do not distinguish between dimers either, counting only their 
number (for example, as showed by Calo [45], the presence of 5-7% dimers leads to the 
appearance of a cluster cascade). 

Thus, we discover new relations, namely an extension of the program "partitio nu- 
merorum" in number theory, from the point of view of the notion of Hartley entropy, 
and indicate possible generalizations of quantization, which lead to an extension of the 
Heisenberg indeterminacy principle [25] . 

Ultrasecond quantization leads to thermodynamics in nanocapillaries and enabled us 
to obtain the superfluidity of liquids in nanotubes [46]- [48], which was confirmed in ex- 
periments (see [49]-[51~]). 

7 Negative pressure and a new critical point of pos- 
sible transition from liquid to "foam" 

First, we use the Wiener quantization of thermodynamics to explain why, in Sec. 5, the 
model of incompressible liquid was used to derive the condition 5*1^=0 = const, where the 
constant is independent of the temperature. 

Because of the Bachinskii relation on the Zeno line, the incompressible liquid model 
leads to a rigid relation between the density (concentration) p and the temperature. 
Since the value of N is undefined for p = (i.e., the concentration p is undefined), the 
temperature is also undefined, and hence, according to the indeterminancy principle, the 
entropy takes a constant value and can be determined. (In addition, we note that, for p = 
0, the activity is equal to 1 at any (undefined) temperature). This precisely means that 
S\n=o = const. As we shall see below, the value of this constant is uniquely determined 
by the new critical point at which the liquid passes into "foam" (dispersion state). The 
value of this point also determines the constant A in the definition of Q, and hence, after 
calculating {p^{V/V CT ), it is possible to obtain the complete description of the distribution 
function for the gas branch in thermodynamics. We note that the indeterminancy at the 
spinodal point agrees with the results of experiments on the absolute instability of the 
spinodal point, which were performed by academician Skripov and his school. 

Now we pass to negative Z. 

As is known, the Bose-Einstein distribution is obtained as the sum of terms of an 
infinitely decreasing geometric progression. If the progression is bounded by the number 



N, then the potential is of the form 

1 -exp^iV 



What is the relationship 17 ! between Ei and 7? 
1) Ultrarelativistic case. E = cp. Here 



E i+ i -Et = cpp 2 dp=\ (p\E i+l ) - p\Ei)^ ~ ~cp 3 (Ei+ l ) = |(^ +1 ) 3 . (65) 



2 

2) Nonrelativistic case. E — f—. Here 



= 1 avs^d 5 _ (^f,r\ a collst(£ ,+ ip . (66) 

3) Consideration of the degrees of freedom. = J^-^- Here 

— Ei = f dp = const (£*fi)(4+^)/(2-Hr). ( 67 ) 
By (ED, 7 = (1 - <r)/(2 + a), and thus 7 < for a > 1. 

As was proved in [TTj [T8j , passing to the limit in the Euler-Maclaurin formula, we 
obtain 

1 r° f 1 iV 



^= (7 + 1)^ + 1) /. i?^T-^T}^ (68) 

In particular, for 7 = —1/2, 

1 /*°° r 1 iv 1 

The absolute value of the derivative of the integrand can readily be estimated by using 
the identities presented below. By the Euler-Maclaurin bounds, this shows that one can 
pass from the sums of the form (1571) to the corresponding integrals with the accuracy 
needed here. 

Hence, denoting N cr = fc , we obtain the following formula for the integral at fi = 0: 



ar( 7 + 2) J e^-l V+ a J e" - 1 
where a = 7 + 1. This implies that 

1 / 1 r 



£i/(i+a) \^ a r(7 + 2) Jo e*- 1 



(71) 



17 The physicists who are not interested in Euler-Maclaurin-type bounds used to pass from sums to 
integrals can skip the following scheme of the proof of these bounds. 
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We obtain 

1 fc 1 . . 1 



e bi _ i e hM -1} d ^ b<* J \et-l £ J ^ 

b«J U e(l + (W2)0/ ba Jo Wl + (W2)OJ 5 

(72) 



Write 



After the change /cq£ = r], we obtain 



^0 / I ^0 ^0 I j tot ^0 



en-l 77(1+77/2)] 6 a y \e"-l 77(1 + 77/2)/ 



en - 1 7/7 7 2(1 + 1) J 6" 6' 

Since rKi+V?) = J ~ 2(1+^/2) > after denotin g ci = / °° 2(1% > we can write 



1 — a «oo j„r> / 7_ \ 1 — a 



, v e ^(i + foy 2 y i + fe V2 

(74) 

Hence, 



b a b a \2 J b a J [ei-l 77(1 - f) J ' 2j 1 + f 6 a 

= -±c+^-c. (75) 
b a b a 

Since &o is the number of particles, b = 1/T, and a = 1 + 7, it follows that /co6 Q 
for 7 > is the value of the Riemann zeta function, £(1 + 7). Therefore, fc 7+1 increases 
for 7 < 1, and the first term of the right-hand side of equation (1751) can be neglected. 
Introducing the function 

M{ 1 + 1) = (c( 7 )/r( 7 + l)) 1/(1+7) , (76) 



we see that the compressibility factor Z 7 = — £(7 + 2)//Vf(7 + l) is subjected to a flexioio 
from 7 > to 7 < 0. 

In this case, we obtain another critical point, which fully corresponds to the physical 
meaning (see [53]). 

Thus, if the compressibility factor is negative, then we divide E/N by T 7+1 with 7 < 
rather than by T, because 



N 



C(7 + 2)7%T = +1 C(7 + 2) 
/W( 7 + l)T red rcd /W( 7 + l) 



18 Since the "Young moduli" for the compression and extension are distinct, a flexion of the spinodal 
occurs. 
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i.e., the energy evaluated for a single particle at fi = (at the "degeneration" point), 
for P < 0, is proportional to Tj^ , i.e., to the temperature taken to a power with an 
exponent less than one. For P < 0, the compressibility factor becomes a dimensional 
quantity; however, this is always considered in this very way on curves in the space{Z, P} 
when using the van der Waals normalization (26). 
We write 

'1 1 



F(0 

For 7 < and /i = 0we have 



j'=i i=i J i=i 

Since the function f(x) = x 1 F(bx) decreases monotonically, we have 

00 00 /»oo /*oo /»oo 

rf J F(6 J ) = V/(j)< / f(x)dx= x^F(bx)dx = b-" / - 1 x~<F{x)dx (78) 

j=l j=l -^0 io JO 

(the Nazaikinskii inequality). 

Thus, the spinodal point j(T) is determined by the relation 

N = b- 1 ((l- 1 ) + 0(b- 1 -i), 7 <0, b=L 

r 

where C is the Riemann zeta function. 

For the critical isotherm we have T r = 1. Hence j(T) = —70- For the point T r and 
the isochor T r 7o+1 ^(7 + 1) we obtain the relation for 7(T r ) from the equation 

T r 70+1 C(7o + 1) = T r C(l + | 7 (T r )|). (79) 

The less T r , the greater |7(T)|. But j(T) cannot exceed 1, because the pressure and Z 
become —00 at this point. 

The negative pressure at the point of the spinodal of liquid phase (for fx — 0) is equal to 
— T 2+7 ^ T ^(2 + 7(T r )), where the value of 7(T) obtained above is negative. The condition 
7(T) < 1 bounds the transition region for the temperature 

T 7o < £( 2 ) — = ^ . (80) 

r "C(7o + l) 6C(7o + l) 1 ] 

For a model we can consider thick rubber with many small slots (see Appendix in [33), 
which is almost incompressible in compression, and because of the holes-slots, has a small 
Young modulus in tension. It is natural to consider tension of a liquid as the appearance 
of nanoholds, i.e., of negative energy PV, which can be treated as negative pressure 
(pressure of the holes). 

Let us now present a graph Fig. 9 for negative pressure for the Lennard- Jones poten- 
tial, where the new critical point is obtained by using a computer experiment!^!. 



19 The absolute zero of temperature is inaccessible. This is visually seen in the logarithmic scale of 
temperatures logT re< i, where the absolute zero corresponds to —00. 




-3 



-4 - 



Figure 9: The spinodal in the coordinates given by the temperature T re( j and the negative 
pressure P rc d- 

8 On Homogeneous Mixtures of Gases 

When considering a gas mixture, we would like to attract the attention at the following 
fundamental point. As is well known, in statistical physics and thermodynamics, the en- 
ergy is sometimes connected with the number of degrees of freedom and the temperature; 
for example, this is the case in the equidistribution law. It turns out in this case that the 
energy depends on the temperature and on the number of degrees of freedom and does 
not depend on the mass. The sequential usage of this conception gave us a continuous 
parameter 7 related to the fractal dimension in the momentum space. Continuing the use 
of this conception in the case of a mixture of pure gases, we should speak of concentration 
rather than on density, i.e., we should neglect the masses of miscible pure gases. 

The sequential application of number theory in thermodynamics, i.e., the consideration 
of the main axiom for the gas mixture (which was in fact made by experimenters, at least 
in the case of air (see [55])), leads to the formulas presented below. 

If two values E\ and E 2 expanded into sums of Aq and N 2 summands, respectively, 
correspond to the fractional dimensions 71 and 72, respectively, and the values of the 
pairs {Si, N x } and {£2,^2} are on the "verge of degeneration," i.e., adding an excessive 
number to Aq and to N 2 leads to the "appearance of the Bose condensate," then, for the 
sum Si + £2 and for N\+ N 2 , any adding of an excessive number to Ni + N 2 also leads to 
the "appearance of the Bose condensate." 

Let pf and p c 2 be the critical concentrations ((in the units cm -3 )), and let Aq and N 2 
be proportional to the molar concentrations, 

Ni/{Ni + N 2 ) = a, (81) 
N 2 /{N X + N 2 ) = (3, (82) 
N = N l + N 2 , a + (3 = 1. (83) 

Since 

E CI = £\, =0 = AVo^ cr w 7 tol TCT = ^ cr (7 cr + 1) 

1 (7 + ^) 

for 7 > 0, it follows that, dividing the equation 

£ s ;=o = 4=o + 4=o = ^ Cr (7i + l)Z?T? + iV 2 cr ( 72 + 1)Z?T? 
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by N™ m , we obtain the relation 

(7 s c u m + l)^su m T s c u r m = a( 7l CT + 1)Z?T? + /3( 72 cr + l)Z^\ (84) 

where Z cr = £(7 + 2)/C t { r ) + 1), the subscripts 1 and 2 refer to the first and second gas of 
the mixture, and the subscript sum refers to the gas mixture. Similarly, 

3 = N(Z-fi/T), 

and, using the additivity of the entropy and dividing by iV^l™, we obtain the relation 

ZZUlZm + 2) = «( 7l cr + 2)Z? + /3( 72 cr + 2)Z?. (85) 

It follows from the above two relations that the quantity 7 = 7 sum almost linearly 
depends also on the values a and T cr = T c s r um . This fact is well known as "Kay's rule" 
in the phenomenological theory of mixtures. For air, we have T c s r um = 232 K, whereas 
T cr = 255 K for oxygen (20% in air) and T cr = 226 K for nitrogen (80% in air). The 
value of T c s r um coincides with the value of this quantity evaluated according to the above 
formulas up to the accuracy of 0.5%. 

He have defined Z CT for ideal gases, and hence also 7 cr for gas mixtures. For a mixture 
of real gases, we must define the function ip™ lx . 

It turns out that, for such mixtures, the Zeno-line is not a segment of a straight line, 
which is the case and is observed experimentally for pure gases. Therefore, for a mixture, 
it is not sufficient to find the values T™ x and p^ lx . One must also define the function 
99™?. This can be carried out by using the following formulas. 

We are interested only in the values of \i\ and // 2 that correspond to the Zeno-line of 
each of the gases (see (l4"])- (j3"8"]) ). 

Since, by assumption, the critical point belongs to the domain of homogeneity, it 
follows that the concentrations a and f3 are preserved, and hence, using equations (61)- 
(63), we obtain the following equation for the sum of the entropies, where k = [ijT: 

(7sum + 2)Z 7sum+ 2(e Ksum ) — K sum 

= «{( 7 i + 2)Z 7l+2 (e Kl ) - Kl } + /3{( 72 + 2)Z 72+2 (e K2 ) - k 2 }, (86) 

where Z 1+2 is equal to the ratio Li 7+2 (e K )/ Li 7+ i(e K ). 

Recall that the values of K\ and k 2 are taken according to the Zeno-lines of the first and 
the second gas, respectively. Hence, using the given values of 7, 71, and 7 2 obtained from 
( IHTT) . we find the value k = n sum , which defines the function (fi^V) by the relation ( l38l . 
This enables us to define the dependence of T c s r um for a mixture of real gases. 

9 Scheme for the evaluation of the percentage for the 
concentration of a mixture of two gases in a liquid 

For the liquid phase, we consider the quasi-ideal case of incompressible liquid. Then, for 
the sum of entropies, we obtain the relation 

T^ +1 { Li 7+2 (e K )( 7 + 2) - «:Li 7+1 (e K )} = T? 1+1 { Li 7l+2 (e Kl )( 7 i + 2) - Li 7l+1 (e Kl )} + 

+T? +1 { Li 72+2 (e K2 )( 72 + 2) - k 2 Li 72+1 (e K2 )}. 

(87) 

20 The full energy does not depend on the masses of the particles, as well as in the theory of Brownian 
particles (see [38]). The latter theory can be expressed in terms of the Wiener path integral, which is 
equal to the Feynman path integral with imaginary Planck constant. This is an additional argument in 
favor of the Wiener quantization. 
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We assume that the state is an equilibrium, i.e., the temperature of the elements of 
the mixture is equal. Then we have 



Ti = T 2 = T. 

If we assume that the full NT-internal energy £ of each of the components is preserved, 
then we obtain 

8 = ^+f>(-^> + T~ { Li7+2 (e«) - k Li^)} 

= | + T' +1 { Li 7+2 (e") - K Li, +1 (e")}. (88) 

The derivative with respect to k of the second term in the expression (|84|) is equal to 
the derivative of Li 7+1 (e K ) with respect to k. 

This means that, integrating over all 7 and taking the variation of (1531) . we obtain the 
condition of chemical equilibrium, 

li5N = firfNi + fi 2 5N 2 . (89) 

This implies a relationship between the chemical potentials, as in the law of mass action. 
A similar formula holds for a mixture of a greater number of components. 

Let us now consider a small correlation sphere for the liquid state of a homogeneous 
mixture of pure substances (as a rule, these are 6-8 molecules) and consider all possible 
replacements of molecules of the gas 1 and the gas 2 in this sphere. Let us form an analog 
of stoichiometric coefficients, as if these were "chemical bonds." 

Now using the relations for chemical potentials, we obtain a system of necessarily 
many equations for the original molecules, as it happens for the law of mass action for 
the original atoms. This enables us to define the mean concentration of each of these two 
gases in the liquid phase. 

Remark 6. Note that it is not rigorously proved in the cycle of papers of the author that 
to any "pure" gas there corresponds a Zeno-line. As is known, for example, the ideal line 
is substantially curved for water under low densities. Water (the creation of dew), as well 
as mercury, must be excluded from our theory of real gases and their mixtures. We have 
heuristically obtained the Zeno-line only for the Lennard- Jones interaction potential. The 
same heuristic proof can be carried our for other interaction potentials. Therefore, one 
can agree that, for pure gases interacting according to the same potential, the Zeno-line 
is approximately a line segment (and this segment is not quite straight in the problem in 
Sec. 1). 

The author thanks D. S. Minenkov for the help in constructing graph 6 and R. V. Nekrasov 
for the help in constructing graph 9. The author is also grateful to Prof. V. S. Borob'ev 
for useful discussions. 

References 

[1] M. A. Leontovich, Introduction to Thermodynamics (GITTL, Moscow-Leningrad, 
1950) [in Russian]. 

[2] H. Poincare, La Science et VHypothe.se (Flammarion, Paris, 1903) [On Science (in 
Russian) (Nauka, Moscow, 1983)] 

39 



[3] H. Temperley, "Statistical mechanics and the partition of numbers: I. The transition 
of liquid helium," Proceedings of the Royal Society of London. Series A, 199 (1058), 
361-375 (1949). 

[4] B. Widom, "Surface Tension and Molecular Correlations Near the Critical Point." 
Journ. Chem. Phys., 43 (11), 3892-3897 (1965). 

[5] L. P. Kadanoff, "Scaling laws for Ising Models Near T c ," Physics, 2, 263-272 (1966). 

[6] K. G. Wilson, "Feynman-Graph Extension for Critical Exponents," Phys. Rev. Lett., 
28 (9), 548-551 (1972).. 

[7] D. Yu. Ivanov, "Critical phenomena in pure liquids," Vestnik SIBGUTI 3, 94-104 
(2009) [in Russian]. 

[8] D. Yu. Ivanov, Critical Behavior of Non-ideal Systems (Wiley- VCH, 2008). 

[9] L. D. Landau and E. M. Lifshits, Theoretical Physics, Vol. 5: Statistical Physics 
(Nauka, Moscow, 1964; Fizmatlit, Moscow, 2003) [in Russian]. 

[10] Yu. G. Pavlenko, Lectures on Theoretical Mechanics (Fizmatlit, Moscow, 2002) [in 
Russian] . 

[11] B. D. Summ, Fundamentals of Colloid Chemistry (Izd. Tsentr "Akademiya," Moscow, 
2007) [in Russian]. 

[12] E. A. Guggenheim, "The Principle of Corresponding States," J. Chem. Phys. 13, 
253-261 (1945). 

[13] V. P. Maslov, "Fluid thermodynamics, an energy redistribution law, a two- 
dimensional condensate, and the T-mapping," Teoret. Mat. Fiz. 161 (3), 420-458 
(2009) [Theoret. Math. Phys. 161 (3), 1681-1713 (2009)]. 

[14] V. P. Maslov, "On Refinement of Several Physical Notions and Solution of the Prob- 
lem of Fluids for Supercritical States", larXiv: 0912.5011v2 [cond-mat.stat-mech]. 

[15] E. M. Apfelbaum, V. S. Vorob'ev, "The confirmation of the critical point-Zeno-line 
similarity set from the numerical modeling data for different interatomic potentials" , 
J. Chem. Phys. 130, 214111, 1-10 (2009). 

[16] A. M. Vershik, "Statistical mechanics of combinatorial partitions, and their limit 
shapes," Funktsional. Anal, i Prilozhen. 30 (2), 19-39 (1996) [Functional Anal. Appl. 
30 (2), 90-105 (1996)]. 

[17] V. P. Maslov and V. E. Nazaikinskii, "On the distribution of integer random variables 
related by a certain linear inequality: I," Mat. Zametki 83 (2), 232-263 (2008) [Math. 
Notes 83 (2), 211-237 (2008)]. 

[18] V. P. Maslov and V. E. Nazaikinskii, "On the distribution of integer random variables 
related by a certain linear inequality: II," Mat. Zametki 83 (3), 381-401 (2008) [Math. 
Notes 83 (3), 345-363 (2008)]. 

[19] E. M. Apfelbaum and V. S. Vorob'ev, "Correspindence between of the ideal Bose 
gas in a space of fractional dimension and a dense nonideal gas according to Maslov 
scheme", Russ. J. Math. Phys. 18 (1), 19-25 (2011). 



[20] L. D. Landau and E. M. Lifshits, Quantum Mechanics (Nauka, Moscow, 1976) [in 
Russian] . 

[21] I. Fisher, The Purchasing Power of Money: Its Determination and Relation to Credit, 
Interest and Crises, (Izd. Delo, Moscow, 2001). 

[22] V. P. Maslov, Zeno-line, Binodal, T — p Diagram and Clusters as a new 
Bose- Condensate Bases on New Global Distributions in Number Theory, arXiv 
1007.4182vl [math-ph] , 23 July 2010. 

[23] V. P. Maslov, "A new approach to probability theory and thermodynamics" , Math. 
Notes, 90, (1), 125-135 (2011). 

[24] V. P. Maslov, "Comparison of the supercritical states of fluids for imperfect gases 
and for a fractal ideal gas," Math. Notes, 87 (3), 303-310 (2010). 

[25] V. P. Maslov, "Tunnel Quantization of Thermodynamics and Critical Exponents", 
Math. Notes, 90, (4), 533-547 (2011). 

[26] V. P. Maslov, "Analytic continuation of asymptotic formulas and axiomatics of ther- 
modynamics and quasithermodynamics," Funktsional. Anal. Prilozhen. 28 (4) 28-41 
(1994). 

[27] V. P. Maslov, "Geometric "quantization" of thermodynamics, and statistical correc- 
tions at critical points," Teoret. Mat. Fiz. 101 (3), 433-441 (1994) [Theoret. Math. 
Phys. 101 (3), 1466-1472 (1995)]. 

[28] V. P. Maslov, "Generalization of Gibbs' axiom, geometric classification of phase 
transitions, and the asymptotics at critical points", Dokl. Ross. Akad. Nauk, 340 
(2), 169-171 (1995). 

[29] V. P. Maslov, "On a class of Lagrangian manifolds corresponding to variational 
problems and problems of control theory and thermodynamics," Funktsional. Anal. 
Prilozhen. 32 (2), 89-91 (1998). 

[30] V. P. Maslov, Quantization of Thermodynamics and Ultrasecond Quantization (Inst. 
Kompyuternykh Issledovanii, Moscow, 2001) [in Russian]. 

[31] V. P. Maslov, Perturbation Theory and Asymptotical Methods (Izd. Moskov. Univ., 
Moscow, 1965; Dunod, Paris, 1972) [in Russian and French]. 

[32] R. P. Feynman and A. R. Hibbs, Quantum Mechanics and Path Integrals (McGraw- 
Hill, New York, 1965 (Dover, Mineola, 2010); Mir, Moscow, 1968). 

[33] V. P. Maslov and P. P. Mosolov, Nonlinear Wave Equations Perturbed by Viscous 
Terms (Walter de Gruyter, Berlin, 2000). 

[34] V. P. Maslov, "Nonstandard characteristics in asymptotic problems," Proceedings 
of the XIX International Mathematical Congress (Warsaw, 1983) [Polish Sci. Press, 
Warsaw, 1984]. 

[35] V. P. Maslov, "Nonstandard characteristics in asymptotic problems," Uspekhi Mat. 
Nauk 38 (6), 3-36 (1983) [Russian Math. Surveys 38 (6), 1-42 (1983)]. 

34 



[36] M. S. Green, "Introduction," Proc. Conference on phenomena in the neighborhood of 
critical points, ed. by M. S. Green & J. V. Sengers (NBS Misc. Publ. 273, Washing- 
ton 1966), pp. ix-xi. 

[37] R. Balescu, Equilibrium and Nonequilibrium Statistical Mechanics, A Wiley- 
Interscience Publication (John Wiley and Sons, New York-London-Sydney-Toronto, 
1975; Mir, Moscow, 1978), Vol. 1. 

[38] V. P. Maslov, "Number-Theoretic Internal Energy for a Gas Mixture" , Russ. J. Math. 
Phys. 18 (2), 163-175 (2011). 

[39] B. van-der-Pol and H. Bremmer, Operational Calculus, Based on the Two-Sided 
Laplace Integral (Cambridge, at the University Press, 1950; Izd. Inostr. Lit., Moscow, 
1952). 

[40] V. P. Maslov and O. Yu. Shvedov, The Complex Germ Method in Many-Particle 
Problems and in Quantum Field Theory (Editorial URSS, Moscow, 2000) [in Russian]. 

[41] V. P. Maslov and S. E. Tariverdiev, "Asymptotics of the Kolmogorov-Feller equation 
for a system with a large number of particles," in Itogi Nauki Tekh., Ser. Probab. 
Theor., Mat. Statist, Theor. Kibern. (VINITI, Moscow, 1982), Vol. 19, pp. 85-120 
[in Russian]. 

[42] V. P. Maslov, "Solution of the Gibbs paradox in the framework of classical mechanics 
(statistical physics) and crystallization of the gas 6*60," Mat. Zametki 83 (5), 787-791 
(2008) [Math. Notes 83 (5), 716-722 (2008)]. 

[43] M. Schonberg, "Application of second quantization methods to the classical statis- 
tical mechanics. I," Nuovo Cim. 9 (12), 1139-1182, (1952). 

[44] M. Schonberg, "Application of second quantization methods to the classical statis- 
tical mechanics. II," Nuovo Cim. 10 (4), 419-472, (1953). 

[45] J. M. Calo. "Dimer formation in supersonic water vapor molecular beams," J. Chem. 
Phys. 62 (12), 4904-4910, (1975). 

[46] V. P. Maslov, "On the superfluidity of classical liquid in nanotubes. I. Case of even 
number of neutrons," Russian J. Math. Phys. 14 (3), 304-318 (2007). 

[47] V. P. Maslov, "On the superfluidity of classical liquid in nanotubes. II. Case of odd 
number of neutrons," Russian J. Math. Phys. 14 (4), 401-412 (2007). 

[48] V. P. Maslov, "On the Superfluidity of Classical Liquid in Nanotubes. Ill," Russian 
J. Math. Phys. 15 (1), 61-65 (2008). 

[49] G. Hummer, J. Rasaiah, and J. Noworyta, "Water conduction through the hydropho- 
bic channel of a carbon nanotube," Nature 414 (8), 188-190 (2001). 

[50] S. Joseph and N. Aluru, "Why are carbon nanotubes fast transporters of water?" 
Nanoletters 8 (2), 452-458 (2008). 

[51] A. Noy, H. Park, F. Fornasiero, et al., "Nanofluidics in carbon nanotubes," Nanoto- 
day 2 (6), 22-29 (2007). 



35 



[52 
[53; 

[54 
[55 



[56 
[57 

[58 

[59 



[60 

[61 
[62 

[63 

[64 

[65 

[66 
[67 



I. A. Molotkov, "Maslov distribution and formulas for the entropy", Russian J. Math. 
Phys. 17 (4), 476-485 (2010). 

K. I. Shmulovich and L. Mercury, "Geochemical Phenomena at Negative Pressures," 
Electronic Scientific Information Journal "Herald of the Department of Earth Sci- 
ences RAS," 1 (24), 1-3 (2006). 

V. P. Maslov, "Theory of chaos and its application to the crisis of debts and the 
origin of the inflation," Russian J. Math. Phys. 16 (1), 103-120 (2009). 

V. V. Sychev, A. A. Vasserman, A. D. Kozlov, G. A. Spiridonov, and V. A. Tsy- 
marnyi, Thermodynamic Properties of Air (Izd-vo Standartov, Moscow, 1978; Hemi- 
sphere Publishing Corp. (National Standard Reference Data Service of the USSR. 
Volume 6), Washington, DC, 1987). 

A. A. Vlasov, "On the vibrational properties of an electron gas," Zh. Eksper. Teoret. 
Fiz. 8, 291-318 (1938) [Sov. Phys. Usp. 10, 721-733 (1968)]. 

W. Wagner, N. Kurzeja, and B. Pieperbeck, "The thermal behaviour of fluid sub- 
stances in the critical region - experiences from recent (p, p, T) measurements on SFq 
with a multi-cell apparatus", Fluid Phase Equilibria 79, 151-174 (1992). 

N. Kurzeja, Th. Tielkes, and W. Wagner, "The nearly classical behavior of a pure 
fluid on the critical isochore very near the critical point under influence of gravity" , 
Int. J. Thermophysics, 20 (2), 531-561, (1999). 

V. P. Maslov, "Taking into account the interaction between particles in a new theory 
of nucleation, quasiparticles, quantization of vortices, and the two-particle distribu- 
tion function," Mat. Zametki 83 (6), 864-879 (2008) [Math. Notes 83 (5-6), 790-803 
(2008)]. 

I. A. Kvasnikov, Thermodynamics and Statistical Physics: Theory of Equilibrium 
Systems (URSS, Moscow, 2002), Vol. 2 [in Russian]. 



http: / /en. wikipedia.org/wiki/Polylogarithm 



L. Lewin, Poly logarithms and Associated Functions (North-Holland Publishing Co., 
New York-Amsterdam, 1981). 

T. L. Hill, Statistical Mechanics: Principles and Selected Applications (McGraw-Hill 
Book Co., Inc., New York-Toronto-London, 1956; IL, Moscow, 1960). 

S. M Stishov, "The Thermodynamics of Melting of Simple Substances," Uspekhi 
Fiz. Nauk 114, 3-40 (1974) [Sov. Phys. Usp. 18, 625-643 (1975)]. 

A. A. Abrikosov, L. P.Gor'kov, I. E. Dzyaloshinskii, Quantum field theoretical methods 
in statistical physics (2ed., Pergamon, 1965). 

B. D. Summ, Basics of colloid chemistry (Academia, Moscow, 2007) [in Russian]. 

V. P. Maslov, "On the Bose condensate in the two-dimensional case, A-point and the 
Thiess-Landau two-liquid model," Teoret. Mat. Fiz. 159 (1), 20-23 (2009). 



[68] V. P. Maslov, "Hypothetic A-Point for Noble Gases," Russ. J. Math. Phys. 17 (4), 
400-413 (2010). 



3fi 



[69] V. P. Maslov, Complex Markov Chains and the Feynman Path Integral for Nonlinear 
Equations (Nauka, Moscow, 1976) [in Russian]. 

[70] N. N. Bogolyubov, On the Theory of Superfluidity, in Selected Works (Naukova 
Dumka, Kiev, 1970), Vol. 2 [in Russian]. 

[71] V. P. Maslov, "On the dependence of the criterion for superfluidity from the radius of 
the capillary," Teoret. Mat. Fiz. 143 (3), 307-327 (2005) [Theoret. and Math. Phys. 
143 (3), 741-759 (2005)]. 

[72] V. P. Maslov, "Resonance between one-particle (Bogoliubov) and two-particle series 
in a superfluid liquid in a capillary," Russ. J. Math. Phys. 12 (3), 369-379 (2005). 

[73] V. P. Maslov, "On the superfluidity of the classical fluid in a nanotube for even and 
odd numbers of neutrons in a molecule," Teoret. Mat. Fiz. 153 (3), 388-408 (2007) 
[Theoret. and Math. Phys. 153 (3), 1677-1696 (2007)]. 

[74] V. P. Maslov, "Quasi-particles associated with Lagrangian manifolds corresponding 
to semiclassical self-consistent fields. Ill," Russ. J. Math. Phys. 3 (2), 271-276 (1995). 

[75] V. P. Maslov, The Complex WKB Method in Nonlinear Equations, (Nauka, Moscow, 
1977) [in Russian]. 

[76] V. P. Maslov, The Complex WKB Method for Nonlinear Equations I (Birkhauser 
Verlag, Basel -Boston-Berlin, 1994). 

[77] V.P.Maslov, "On the integral equation 

u{x) = F(x) + /G(x,0M+~ 2)/2 (0 d C// u +~ 2)/2 (0 d C forn = 2 and n = 3", Mat 
Zametki, 55 (3), 96-108, (1994); Math. Notes, 55 (3), 302-311, (1994). 

[78] V.P.Maslov, "Spectral series, superfluidity, and high-temperature superconductiv- 
ity", Mat. Zametki, 55 (3), 1349-1352, (1995). 

[79] V.P.Maslov and O.Yu. Shvedov, "On the number of particles in the condensate of 
weakly ideal Boze gas", Mat. Zametki 61 (5) 790-792, (1997). 

[80] V.P.Maslov, "On an averaging method for the quantum many-body problem", Funk- 
tsional. Anal. Prilozhen. 33 (4), 50-61, (1999). 

[81] V. P. Maslov, "Mathematical aspects of weakly nonideal bose and fermi gases on a 
crystal base", Funktsional. Anal, i Prilozhen. 37 (2), 16-27 (2003). 

[82] V. P. Maslov, "Demonstrativeness in Mathematics and Physics,", Russ. J. Math. 
Phys.l9(2), 163-175(2012). 

[83] P. Erdos, "On some asymptotic formulas in the theory of partitions," Bull. Amer. 
Math. Soc. 52, 185-188, (1946). 

[84] P. Erdos and J. Lehner, "The distribution of the number of summands in the parti- 
tions of a positive integer," Duke Mathematical Journal 8 (2), 335-345 (June 1941). 



37 



